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Abstract: Construction duration and schedule robustness are of great importance to ensure efficient 
construction. However, the current literature has neglected the importance of schedule robustness. 
Relatively little attention has been paid to schedule robustness via deviation of an activity’s starting 
time, which does not consider schedule robustness via structural deviation caused by the logical 
relationships among activities. This leads to a possibility of deviation between the planned schedule 
and the actual situation. Thus, an optimization model of construction duration and schedule robustness 
is proposed to solve this problem. Firstly, duration and two robustness criteria including starting time 
deviation and structural deviation were selected as the optimization objectives. Secondly, critical chain 
method and starting time criticality (STC) method were adopted to allocate buffers to the schedule in 
order to generate alternative schedules for optimization. Thirdly, hybrid grey wolf optimizer with 
sine cosine algorithm (HGWOSCA) was proposed to solve the optimization model. The movement 
directions and speed of grey wolf optimizer (GWO) was improved by sine cosine algorithm (SCA) 
so that the algorithm’s performance of convergence, diversity, accuracy, and distribution improved. 
Finally, an underground power station in China was used for a case study, by which the applicability 
and advantages of the proposed model were proved. 


Keywords: optimization model; construction duration; schedule robustness; critical chain method; 
STC method; hybrid grey wolf optimizer with sine cosine algorithm 


1. Introduction 


Effective schedule management is becoming increasingly significant because of its social and 
economic benefits. Optimization is an important issue in schedule management. Both construction 
duration and schedule robustness, the latter of which is the schedule’s ability to resist duration 
increases due to interference, are necessary to be taken as the optimization objectives. However, 
schedule robustness is rarely considered as an optimization objective in current research, which causes 
deviations between the actual schedule and the planned one. Therefore, it is necessary to conduct 
optimization of construction duration and schedule robustness. 

The basics of optimization is to determine objectives. Besides the objective of construction duration, 
it is also necessary to determine appropriate robustness criteria as optimization objectives. A few 
general robustness criteria have been defined to measure schedule robustness, including the summation 
of each activity’s free float [1], ratio of each activity’s free float and duration [2,3], linear weighted 
summation of the total free float [4], and linear weighted summation of the time deviations [5]. 
For complex construction projects, the schedule is established by breaking down the construction 
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project into basic units (which are called activities) and determining the starting time of activities and 
logical relationships among them. Starting time of activities and logical relationships among them are 
two indispensable parts, thus the schedule robustness should be measured from these two aspects. 
But, in the very few multi-objective optimization studies considering both construction duration and 
schedule robustness, only the starting time deviation is taken as an objective [6]. The neglect of schedule 
robustness from the perspective of logical relationships among activities—that is, when structural 
deviation isn’t taken as an objective—will cause the construction site to be chaotic. Therefore, it is 
practical and necessary to take construction duration, starting time deviation, and structural deviation 
as the optimization objectives. 

The method to improve schedule robustness is a necessary condition to generate alternative 
schedules for optimization. Buffer allocation methods are proven effective methods to improve schedule 
robustness. These methods are mainly based on two ideas: Allocating buffers centrally at the end of 
the schedule chain and allocating buffers scatteredly to the schedule. The critical chain method [7] 
is the main method to allocate buffers centrally [8-10] and its advantages have been verified [11-13]. 
The methods of allocating buffers scatteredly mainly include virtual activity duration extension 
(VADE) [14], resource flow dependent float factor (RFDFF) [15] and starting time criticality (STC) [16], 
which are also widely studied and applied. These aforementioned studies indicate that allocating 
buffers centrally improves the robustness of the whole project and allocating buffers scatteredly 
improves the robustness of activities within the project. Researchers tend to allocate buffers centrally or 
scatteredly according to the research emphasis [17], rather than allocating both kinds of buffers at the 
same time. However, for complex construction projects, the orderly construction of each activity and 
the timely completion of the whole project are equally important. Therefore, it is preferable to allocate 
both kinds of buffers to the schedule in proper proportion, so that the schedules for optimization of 
construction duration and schedule robustness can be generated. 

Metaheuristic algorithm is an effective tool for solving optimization problems. Several 
metaheuristic algorithms have been developed, such as particle swarm optimization (PSO) [18], 
evolutionary algorithm (EA) [19], multi-verse optimization (MVO) [20], genetic algorithm (GA) [21], ant 
colony optimization (ACO) [22], and whale optimization algorithm (WOA) [23]. GWOis a very efficient 
metaheuristic algorithm due to its simple mathematical model [24]. In GWO, the search direction is 
determined by the best three wolves and followed by other candidate wolves. It has been successfully 
applied for solving optimization problems in the fields of flow shop scheduling [25-27], computer 
science [28-30], mathematics [31-33], water resources [34-36], energy [37-39], and so on. But this kind of 
exploitation mechanism makes it prone to stagnation in a local optimum in multi-objective optimization. 
The sine cosine algorithm (SCA) [40], which has been widely applied in the fields of engineering [41-43], 
computer science [44—46], control system [47,48], energy [49-51], and instrument [52,53], requires 
the solutions to fluctuate towards or outwards of the true Pareto optimal front based on the sine 
cosine function. SCA is superior in balancing exploration and exploitation, which makes it able to 
approximate the true Pareto optimal front. Therefore, SCA is useful to improve the searching capability 
of multi-objective GWO so that the global optimal solutions of the optimization model of construction 
duration and schedule robustness can be obtained. 

In summary, the current studies have paid little attention to schedule robustness via starting time 
deviation, which neglected the schedule robustness via structural deviation. Additionally, a more 
efficient algorithm is worth being developed to solve the optimization model. Therefore, the aim of this 
research is to propose a novel optimization model of construction duration and schedule robustness 
based on a hybrid grey wolf optimizer with the sine cosine algorithm (HGWOSCA), which can be 
divided into three detailed points: 


1. Construction duration, starting time deviation, and structural deviation are taken as the objectives 
of the optimization model. 

2. Buffers are centrally allocated to the schedule using the critical chain method and scatteredly 
using the STC method so as to generate alternative schedules for optimization. 
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3. The hybrid grey wolf optimizer with sine cosine algorithm (HGWOSCA) is proposed to solve the 
optimization model. 


The remainder of this paper is organized as follows: The framework of this paper is shown in 
Section 2. Section 3 describes the details of the proposed methodology. Section 4 presents a case study 
and results. Section 5 provides the discussion to imply the superiority of the proposed optimization 
model. Section 6 highlights the conclusions of this paper. 


2. Research Framework 


In this research, an optimization model of construction duration and schedule robustness is 
proposed so that a schedule with both short construction duration and good robustness can be 
obtained. The research framework is shown in Figure 1. Starting time deviation rs, structural deviation 
rp, and construction duration T were the optimization objectives. Critical chain method and STC 
method were adopted to allocate project buffer (PB), feeding buffer (FB), and scattered buffers (SB) 
to the schedule so as to generate the schedules for optimization (see Appendix A). Hybrid grey wolf 
optimizer with sine cosine algorithm (HGWOSCA), where the position updating way of WGO is 
improved by SCA, was proposed to solve the optimization model. Then, an underground power 
station in China was used as a case study and the applicability and advantages of the proposed 
optimization model are discussed and verified. 


Optimization model 


* Optimized schedules 


max|s,,+T(i-1)]<S,,i21 No. PB/day FB/day  SB/day r/Jday r Tiday 

w, (|P; - p,|) Subject to: į max[ s, , +T(i-1)+SB+FB]<5,,i21 75.12 0.6027 2398 

‘ 76.13 0.6072 2393 
75.05 0.6027 2401 

75.28 0.6126 2396 

5 75.07 0.6027 2399 


* Gantt chart of optimized Schedule No. 1 MEE initial schedule 
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Figure 1. Research framework. 


3. Methodology 


This section presents the description of the optimization model of construction duration and 
schedule robustness. In addition, the methodologies of robustness criteria, generation of alternative 
schedules for optimization, and the proposed HGWOSCA are introduced. 
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3.1. Problem Formulation 


The optimization model of construction duration and schedule robustness is expressed as follows: 


Opt (Ts, fp, T) 
n 
fs = 2 w;(|S; — sil) 
1= 


r 1 
lp = X w;(|P; — pil) ( ) 
{= 
T= > (SB + T(i)) + PB 
1€critical path 
subject to: 
max(|s;_1+T(i-1)] < S,i2 1 
max|s;_7+T(i-1) + SB + FB] < S;,i 21 (2) 


yir<R 


where, three optimization objectives are proposed. rs is the robustness criterion proposed from the 
aspect of deviation of an activity’s starting time, rp is the robustness criterion proposed from the 
aspect of logical relationships among activities, and T represents the total construction duration of the 
schedule. The detailed descriptions of the symbols in this paper are shown in Appendix A. 


3.2. Robustness Criteria 


Robustness criteria are the foundation for the optimization model. Since activities and logical 
relationships among activities are two indispensable components of a construction schedule [6], 
two criteria (i.e., starting time deviation (rs) and structural deviation (rp)), are proposed to measure 
the schedule robustness from these two aspects. r, is the robustness criterion proposed from the 
aspect of activities’ starting time, and ry is the robustness criterion proposed from the aspect of logical 
relationships among activities. These two criteria provide a more comprehensive understanding of 
measuring schedule robustness. 

The starting time deviation r; is the absolute value of the weighted summation of the deviations 
between the starting times of the actual activities and the planned activities. 


rs =) wil|Si- si). (3) 
i=0 


The structural deviation ry is the absolute value of the weighted summation of the deviations 
between the orders of the actual activity and the planned activity. 


tp =) w;(|Pi-pil). (4) 
i=0 


As seen from the above equations, the smaller the r; and ry values, the better the robustness. 


3.3. Generation of Alternative Schedules for Optimization 


With the robustness criteria as the optimization objectives, the next step was to improve schedule 
robustness to generate alternative schedules for optimization. This included two parts (as shown in 
Figure 2): Identifying the critical chain in the schedule by critical chain method and allocating buffers 
to it by using the STC method. 
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Figure 2. Generation of alternative schedules with buffers allocated based on critical chain method and 
the starting time criticality (STC) method. 


The critical chain is a chain composed of activities that are interdependent due to time conflicts 
and resource conflicts. Since the constraints which exist in each schedule correspond to the bottleneck 
elements that prevent the project from being completed on schedule [7], it is effective to identify the 
conflicting activities so as to improve the schedule robustness. 

The identification process of a critical chain was as follows: 


1. Compute the initial schedule via simulation. 

2. Identify the activities that can be postponed on the condition that the schedule fulfils the resource 
constraints and logical relationships. An activity which cannot be postponed was denoted 
as NonAi; otherwise, it was denoted as MovAi. A MovAi which does not satisfy the resource 
constraints was denoted as NonBi; otherwise, it was denoted as MovB1. 

3. The critical chain was identified by combining NonAi and NonBi. If there was more than one 
critical chain, the chain with the most activities was selected as the final critical chain, and the 
others were treated as non-critical chains. 


After the critical chain was identified, three types of buffers, including project buffer (PB), scattered 
buffer (SB), and feeding buffer (FB), were allocated to it. PB was allocated at the end of the critical chain 
to ensure the project was completed on time. SBs were allocated in front of highly critical activities to 
ensure that these activities started on schedule. FBs were allocated where non-critical chains feed into 
the critical chain in order to prevent a delay of the activity on the critical chain when the activities on 
the non-critical chain were delayed. 

The sizes of FB and the sum of PB and SB were determined using the elasticity coefficient method 
which was improved by the Dezert-Smarandache Theory (DSmT) [54,55]. 

In the elasticity coefficient method, the (a, b, c) in the three-point estimation method is replaced by 
(a, b) and c is considered in the elasticity coefficient: K = (c —a)/(b —a). In this way, the impacts of 
interference factors on schedule can be reflected in buffers. The closer c approaches a, the less likely it 
is to delay the activity and the smaller the K value. The closer c approaches b, the more likely it is to 
delay the activity and the greater the K value. Therefore, the buffers can be calculated through K. 


FB = Y \{(b¢ -a°) x K), (5) 
ceB 


where B denotes the set of activities related to the buffers. 

Additionally, in the elasticity coefficient method, a is the activity duration without considering the 
impacts of interference factors and is calculated by simulation. The activity durations are represented 
by b and c and consider varying degrees of impacts of interference factors and are calculated as: 


b =a x (1+min(m;(01), m;(02), m;(03))); (6) 


3 
c=ax (1+ Y mj(8)), (7) 
n=1 
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where m;(@;,) is computed using the DSmT method. Delay, downtime, and breakdown happen 
during construction as a result of interference factors. Hence, it is reasonable to determine the 
size of buffers according to the impact of interference on schedule. DSmT recognizes the fusion of 
occurrence probabilities of interference factors. It reasonably solves the problem of large conflicts 
among interference factors and achieves the unification of occurrence probabilities and impacts of 
interference factors on schedule. 

Interference factor j (j = 1 — 6) represents the geological factor, technical factor, environmental 
factor, management factor, material factor, and the mechanical factor, respectively. The construction 
parameter k (k = 1-11) represents surveying setting out, charge explosion, safe disposal, explosion slag 
cleaning, ventilation, disposal of unfavorable geology, other downtime, hand drill, three arm trolley, 
dump truck and loader, and the consequence of interference factors, respectively. Factor h (h = 1-3), 
represents delay, downtime, and breakdown, respectively. Factor m;(0;) is defined as: 


m;(01) =m,(y1) ®... 8m, (y6) 
m;(02) =m,(y7)®..8mj(yo) ; (8) 
m;(03) = m,(y10) ®m;(y11) 

VA #2 EDU, [m +m](A) = D mı(01)m2(02). (9) 


01,02€DU,6, NOo=A 


Then, the STC method and golden section rule were used to determine the locations of scattered 
buffers. The process is as follows: 


1. The STC method was adopted to measure the necessity of allocating the buffer before an activity. 
T 
stc(i) = X m;(On) X wi, (10) 
n=1 


3 
where 4 L m;(;,) denotes the possibility that the activity i does not start on schedule. 
n= 


2. Based on the golden section rule, 38.2% of the buffers that need to be allocated on critical chain 
were allocated at the end of the critical chain as the project buffer. Based on this same rule, 61.8% 
of the buffers that need to be allocated on critical chain were allocated in front of activities as the 
scattered buffers. 


PB = 0.382) {(b° — a°) x K). (11) 
ceB 

SB = 0.618) {(bf — a°) x KẸ}. (12) 
ceB 


3. Sort the activities on critical chain from high to low in order of criticality. Add their buffers one 
by one. 

4. Judge whether the sum has reached 61.8% of the buffers. It not, go back to step 3. 

5. Ifthe sum has reached 61.8% of the buffers, the added activities are the activities that scattered 
buffers are allocated before. 


3.4. Optimization of Construction Duration and Schedule Robustness Based on HGWOSCA 


In this section, a hybrid multi-objective grey wolf optimizer with sine cosine algorithm is applied 
to solve the optimization model. The movement directions and speed of GWO was improved by SCA. 
Then the validity and superiority of the HGWOSCA was demonstrated by comparison with four other 
algorithms. The solution flow of HGWOSCA for T-r;-rp optimization model is presented. 
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3.4.1. Theory of HGWOSCA 


The HGWOSCA contains three phases: Encircling, hunting behavior improved by SCA, 
and modified searching and attacking behavior. Four types of grey wolves named alpha (a), beta (£), 
delta (ô), and omega (w) were employed for mathematically modeling their social hierarchy. The first 
three fittest solutions are a, fp, and 6, respectively. The rest of the candidate solutions are w. In GWO, 
the optimization is guided by a, f, and 6, and the w wolves follow these three. 

The encircling behavior is simulated by the following equations: 


D =|C-X,(t) -X(®), (13) 
X(t+1) = X(t) - AD, (14) 


> > 
where t represents the current iteration, X, is the position vector of the prey, and X represents the 


> > 
position vector of a grey wolf, A and C are coefficient vectors, which are calculated as follows: 


A=2a-n-4a, (15) 
C=2-7r, (16) 


where 11, r are random vectors in [0, 1]. 
The hunting behavior is simulated by the following equations: 


i 


TID 
QI QL Ol 
SSL RAL PSI 
| 
xl dl l 


D, (17) 
= E 
X41 — Xa — Á] . (Da) 
Xp = Xg - A2 : (Dg) » (18) 
3 = 446 T 13: ( 5) 
> X +X +X, 
X(t+1) = arara (19) 


GWO was originally developed for solving continuous optimization problems and cannot be 
used to solve the combinatorial optimization problems in this paper. Thus, the search operator was 
modified as follows: 


. shi ft(X(t),Xa-X(t)), if r< 4 
X(E+1) =} shift(X(t),Xp— X(t), ifk<r<2 . (20) 
shi ft(X(t), X; — X(t), otherwise 


If the units were beyond the boundary range, the boundary element was connected to it in the 
other direction. 

The position of a is determined by the position updating equation of SCA [40], by which the 
position, speed, and convergence accuracy of the grey wolf agent are improved and the convergence 
performance of GWO is extended. 





, KA< 0.5 


ae Xi +r x sin (r2) x [r3 P! A 21) 
l , 420.5 


XÍ + r1 X cos (r2) X [r3 Pi =X 
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where, r1, r2, and r3 are random numbers, and r4 is a random number used for exploitation and 
exploration. 


3.4.2. Performance of the Proposed HGWOSCA 


The efficiency of the proposed HGWOSCA was tested via the patterson_patl and patterson_pat20 
in the well-known PSPLIB data set [56], as shown in Figure 3. Multi-objective particle swarm 
optimization (MOPSO), multi-objective grey wolf optimizer (MOGWO), multi-objective evolutionary 
algorithm (MOEA), and multi-objective multi-verse optimization (MOMVO) were selected for 
comparison with the HGWOSCA. The tests were implemented on a computer with Intel Core 
i5, 2.39 GHz, 131,072 MB RAM, with a Windows 10 operating system. 
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Figure 3. Network of patterson_patl and patterson_pat20 [56]. 


Four metrics including generational distance (GD), inverse generational distance (IGD), ratio of 
non-dominated individual (RNI), and Spread were employed. 

(1) GD represents the distance between the Pareto front and the optimal Pareto front and is 
formulated as: 


GD = a (22) 


where PF is the number of solutions in the Pareto front. The lower the GD value, the better the 
convergence performance of an algorithm. 

(2) IGD represents the distances between each solution in optimal Pareto front and is formulated 
as: 


2 
2 x € OPF dp 


where OPF is the number of solutions in optimal Pareto front. The lower the IGD value, the better the 
convergence performance of an algorithm. 
(3) RNI represents the ratio of the number of non-dominated individuals to the size of the 


population and is formulated as: 
Nuon 


RNI = N x 100%, (24) 
where Nņpnon is the number of non-dominated individuals and N is the number of all individuals. 
The higher the RNI value, the better the solution quality. 

(4) Spread represents the extent of spread achieved among the front and is formulated as: 


PF — 
3 han +E ka a-a 
Spread = ————————, (25) 


no g J 
y Mar dé + |PF| -d 
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where ng is the number of objectives. The lower the spread value, the better the distribution of solutions. 

Figures 4 and 5 show the four metrics on the algorithms mentioned above over 50 runs. 
The GD value, IGD value, and spread value of HGWOSCA were the lowest among the five 
algorithms, and the RNI value of HGWOSCA was the highest among the five algorithms. This result 
indicates the good performance of the proposed HGWOSCA. Thus, the HGWOSCA is a good fit for 
multi-objective optimization. 
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3.4.3. Solution Flow of HGWOSCA for T-rs-rp Optimization Model 


The proposed HGWOSCA was applied to solve the optimization model of construction duration 
T, starting time deviation rs, and structural deviation rp. A basic flowchart of this is shown in Figure 6. 


Encoding and decoding for HGWOSCA 


Initialization: 
a. One initial schedule is generated using critical chain method; 
b. Other schedules are generated randomly. 


...The search operator is modified to solve 


Search for the prey z optimization of discrete variables 


The position updating way of a is improved by SCA 


Compute the three criteria for each schedule and sort the schedules using Pareto. 
Evaluate newly generated individuals. 


Combine previous population and newly generated 


individuals to form the new population. 











Sort the population according to the social hierarchy 
and select the good population as the next population. 


t=t+1 t< maximum iteration? 


No 


Output the best solutions 


Figure 6. A flowchart of the hybrid grey wolf optimizer with sine cosine algorithm (HGWOSCA). 


Encode the schedule. 
Initialization: To ensure a good diversity of solutions, one schedule was generated using the 
method proposed in Section 3.3 and other schedules were generated randomly. 


NO e 
Ne” N 


(3) Search for the prey. The search operator of GWO was originally developed to solve continuous 
optimization problems [57], which was not suitable to solve the optimization problems of discrete 
variables. Thus, in this paper, the search operator was modified as shown in Section 3.4.1. 
Secondly, the position updating way of a was improved by SCA to improve the searching 
capability of GWO as shown in Section 3.4.1. 
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(4) Compute the three criteria for each schedule and sort the schedules using Pareto relationship. 
According to the Pareto dominance nature of multi-objective optimization problems, there is 
no unique optimal solution but a set of non-dominated solutions. Thus, the population can be 
divided into several ranks according to the Pareto dominance relationship. The solutions at the 
first level rank are denoted as a solutions, the solutions at the second level rank are denoted as £ 
solutions, the solutions at the third level rank are denoted as 6 solutions, and the other solutions 
are donated as w solutions. 

(5) Evaluate: The newly generated individuals were evaluated based on the fitness values. The best 
individuals were saved and combined with the previous population to form the new population. 
Then the populations were sorted according to the social hierarchy and the good population was 
selected as the new population. 

(6) Judgment: Determine whether the algorithm satisfies the termination condition. If yes, output 
the best solutions; otherwise, let t = t+1 and go back to step 5. 

(7) Output: Output the best solutions which are the schedules with both short construction duration 
and good schedule robustness. 


4. Case Study 


A real-life underground power station in Southwest China was selected as a case study to 
demonstrate the rationality of the research in this paper. The layout of the underground power station 
is shown in Figure 7. It is composed of a main workshop, main transformer, last chamber, and several 
branch tunnels. 






ast chamber 


Tailrace tunnel 














Branch tunnel 


Main transformer 


Diversion tunnel 


Main workshop 


























Figure 7. Layout of the project. 


According to the construction arrangement provided by the design institute, the initial schedule’s 
information was: T = 2376 days, rs = 82.66 days, and r, = 0.98. Then, the optimized schedules were 
obtained after 200 iterations of optimization. Five of them are shown in Table 1. These schedules had 
better performance for both construction duration and schedule robustness. The optimized schedules 
corresponded to a set of Pareto optimal solutions rather than a unique solution. This is because a 
conflict exists between construction duration and schedule robustness, and they cannot be directly 
compared. It is impossible to improve either objective without weakening the others. Therefore, the 
optimized schedules in Table 1 met the requirements of multiple objectives and provided a variety of 
satisfying schedules. Gantt chart of main workshop layer #1—#4 of the optimized schedule No. 1 in 
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Table 1 is shown in Figure 8 as a detailed presentation, in which the duration of each activity and the 
schedule of robustness is clearly shown. 


Table 1. Optimized schedule. 


Schedule Number PB/day FB/day SB/day _ r,/day fp T/day 


1 19 22 14 75.12 0.6027 2398 
2 19 22 6 76.13 0.6072 2393 
3 19 22 30 75.05 0.6027 2401 
4 20 22 12 75.28 0.6126 2396 
5 20 22 23 75.07 0.6027 2399 





301: Main workshop layer #1 - diversion tunnel, T = 56 days MSS Initial schedule 

302: Main workshop layer #2 - diversion tunnel, T=55 days |E Optimized schedule 
303: Main workshop reserved protective layer , T = 33 days 
304: Main workshop layer #1 - downstream expanding excavation , T = 42 days 
305: Main workshop lagging , T = 17 days 
306: Main workshop layer #1 - upstream expanding excavation , T = 43 days 
308: Main workshop layer #3-1-1 , T = 33 days 
309: Main workshop layer #3-1-2 upstream , T = 15 days 
310: Main workshop layer #3-1-2 downstream , T = 15 days 
311: Main workshop layer #3-2-1 , T = 34 days 
pmm ggg o aa h 312: Main workshop layer #3-2-2 upstream , T= 16 days 


Ts16-intiat = 48 days 313: Main workshop layer #3-2-2 downstream , T = 16 days 
Tete optimised = 50 days 







314: Main workshop layer #3-3 upstream , T = 14 days 


316: Main workshop T317-imtial = 48 days 315: Main workshop layer #3-3 downstream , T = 14 days 


layer #4-1 1 317-optimized = 49 days k = 316: Main workshop layer #4-1, T=50 days | 
317: Main workshop T318. intial = 66 days iN 317: Main workshop layer #4-2, T=49 days' 
layer #4-2 Saas) T318-optimized = 68 days: 318: Concrete of main ! 





Concrete of main workshop : | workshop, T= 68 days 


Figure 8. Gantt chart of the main workshop layer #1—#4 of the optimized schedule No. 1. 


5. Discussion 


To determine the effectiveness of the proposed HGWOSCA for optimization of construction 
duration and schedule robustness, four algorithms including MOPSO, MOGWO, MOEA, and MOMVO 
were selected for comparison. To make a fair comparison among these algorithms, the same encoding 
scheme and initialization were used as HGWOSCA. Each algorithm ran 500 iterations independently 
and the maximum number of solutions was set as 200. 

The Pareto front obtained by the five algorithms above is shown in Figure 9, in which (a) provides 
a three-dimensional view with rs, rp, and T; (b) provides a two-dimensional view with rs and T; 
(c) provides a two-dimensional view with r, and T; and (d) provides a two-dimensional view with 
rs and rp. It can be seen that the HGWOSCA can approach closer to the Pareto front than the 
other algorithms, which implies that the HGWOSCA can obtain better solutions than the other four 
algorithms. The proposed HGWOSCA is therefore superior to its compared algorithms for optimization 
of construction duration and schedule robustness. 
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Figure 9. Pareto fronts (PF) obtained by different algorithms under different angles: (a) PF with rs, rp 
and T; (b) PF with rs and T; (c) PF with rs and rp; (d) PF with rs and T. 


The four metrics of HGWOSCA and the other four algorithms used for comparison are shown 
in Table 2. The GD value of HGWOSCA was the lowest, which indicates that HGWOSCA had a 
better convergence performance. The IGD value of HGWOSCA was the lowest, which indicates that 
HGWOSCA had the better diversity and convergence performance. The RNI value of HGWOSCA 
was the highest, which indicates that HGWOSCA had the better accuracy performance. The spread 
value of HGWOSCA was the lowest, which indicates that HGWOSCA had the better distribution 
performance. Therefore, HGWOSCA is superior in its performance of convergence, diversity, accuracy, 
and distribution for optimization of construction duration and schedule robustness. 


Table 2. Four metrics based on different algorithms. 


Algorithm 
MOPSO MOGWO MOMVO MOEA HGWOSCA 
Metric 
GD 0.0636 0.0518 0.1599 0.0529 0.0500 
IGD 2.9819 0.8866 16.361 0.9196 0.7176 
RNI 0.2079 0.0528 0.0792 0.0660 0.597 
Spread 0.7186 0.7164 0.9680 0.8096 0.6984 


6. Conclusions 


Construction duration and schedule robustness have great significance in construction schedule 
management. However, most related researches which take both construction duration and schedule 
robustness into account only considered schedule robustness via starting time deviation. Additionally, 
a more efficient algorithm is worth being developed to solve the optimization model. To overcome these 
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limitations, an optimization model of construction duration and schedule robustness was proposed, 
which considered starting time deviation and structural deviation. An underground power station in 
China was used for a case study to verify its applicability and advantages. The major contributions are 
summarized as follows: 


e Considering schedule robustness via starting time deviation and structural deviation, a novel 
optimization model of construction duration and schedule robustness was proposed. 

e Alternative schedules were generated using critical chain method and STC method to allocate 
buffers centrally and scatteredly. 

e The hybrid grey wolf optimizer with sine cosine algorithm (HGWOSCA) was developed to solve 
the optimization model. 


Overall, the proposed optimization model of construction duration and schedule robustness 
is applicable and has practical significance for multi-objective construction schedule management. 
In future research, an optimization system version centered on human-computer interaction can be 
generated to improve solving efficiency and operability for optimization. 
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Appendix A Symbol Descriptions 





Symbol Definition Symbol Definition 
i Robustiices criterion, starting time Probability of change of construction 

deviation EERS a7 miyi) parameter k in the case of interference factor 
fp Robustness criterion, structural deviation j occurs 
T Total construction duration mlO) The impacts of interference factor j on 
1 Activity number, 7 = 1,2,...,n ae schedule 
Wj Denotes the weight of activity í nu The set generated by intersection and union 

eas eee f el i 
Si Actual starting time of activity í PE Cement 1H 
stc(i) Starting time criticality of activity i 
Si Planned starting time of activity 7 
ses The first, second, and third best solution in 
P; Actual order of activity i a, B, Oo . 
ee GWO, respectively 
Di Panned order of activity 7 
T (1) Construction duration of activity 7 w The rest of solutions in GWO 
l _ 3 3 3 OS Vector or variable calculated by vectors in 
NonAi An activity which cannot be postponed Ari TD Xe GWO y 
MovAi An activity which can be postponed xt Position of the current solution in /th 
l An MovAi which does not satisfy the ; dimension at tth iteration 

NonBi wy SA: ae 

resource constraints pt Position of the destination point in Ith 

dimension 

An MovAi which satisfies the resource ; 
MovBi cid i ai i The Euclidean distance between each 

constraints ee 

n 7 j _ dy solution in the pth front and the nearest 

K Elasticity coefficient 


Interference factor number 


member in the optimal Pareto front. 
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k Construction parameter number i an cone 
The Euclidean distance between the gth 
h Consequence number of interference factors 4, extreme solutions and the solutions in the 
Optimistic, pessimistic and the most likely PF 
ane duration of an activity, respectively 
Acronyms 
Symbol Definition Symbol Definition 
PB Project buffer GWO Grey wolf optimizer 
FB Feeding buffer SCA Sine cosine algorithm 
= Ao MOEA Multi-objective evolutionary algorithm 
IGD Inverse generational distance MOCWO Mulzobjective grey wolf opamizer 
RNI Ratio of non-dominated individual 
Spread The extent of spread achieved among the MOPSO Multi-objective particle swarm optimization 
PE ae a MOMVO Multi-objective multi-verse optimization 
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